Copyright © 1976 American Telephone and Telegraph Company 

The Bell System Technical Journal 

Vol. 55, No. 4, April 1976 

Printed in U.S.A. 



Digital Single-Tone Generator-Detectors 

By R. P. KURSHAN and B. GOPINATH 

(Manuscript received December 5, 1975) 

A class of digital, linear generator-detectors, based upon cyclotomic 
polynomials, which have simple implementation and operate without 
roundoff errors, is proposed. It is shown how these filters are optimal among 
all linear generator-detectors which have no roundoff required in the feed- 
back loop. The complexity of various cyclotomic filters are compared. 
These filters in general require far fewer binary adds/s than conventional 
second-order filters used for the same purpose. 

I. INTRODUCTION 

Devices for pure tone generation and detection have widespread 
applications. The most notable examples are Touch-Tone® signaling, 
frequency shift keying (fsk), and multifrequency (mp) signaling. 
Associated with such devices are problems of stability and predict- 
ability, which in practice are dealt with on an individual basis, using 
techniques peculiar to the particular application. When these devices 
are realized digitally, the above problems are manifest from errors due 
to operational roundoff. 

Generally, tones for signaling are analog signals of the form A sin tot 
(A is the amplitude, 2ir/u is the period, and u/2ir is the frequency). 
Devices that generate these tones are usually oscillators of various 
kinds. Because of the requirement of structural stability, in practice 
these devices are limit cycle oscillators. These are simulations and 
realizations in hardware of nonlinear differential equations that have 
limit cycles. Because of the complexity of these equations, the ampli- 
tude and frequency are not easily predicted from given values of resis- 
tors and capacitors in the network. 

For detection of these tones, linear analog filters are frequently used. 
These are also used as generators, when the duration of the signal is not 
too long compared to the period. However, passive linear analog oscil- 
lators require inductors which are bulky, and the frequency and 
amplitude of these oscillators can vary with changes in value of the 
inductors and capacitors due to environmental conditions. Active 

469 



linear oscillators using rc elements are used in many applications. 
However, they also generally need some form of limiting and end up 
being nonlinear devices, thus usually preventing them from being used 
as receivers. 

Digital oscillators, on the other hand, are almost insensitive to chang- 
ing parameter values and produce stable repeatable waveforms. How- 
ever, in the mechanization of these oscillators (which are usually 
based upon second-order linear equations) , roundoff in multiplication 
and addition produce errors in the feedback that lead to limit cycles 
and can significantly impair the signal quality. Also, when such linear 
digital devices are used as receivers, the precision required for satis- 
factory performance goes up quite rapidly with increasing Q. Although 
the effects of this can be satisfactorily controlled in certain specific 
applications (see, for example, Ref. 1), the difficulties, in general, can- 
not be ameliorated except by increasing the accuracy of computations. 2 

In this paper, we present a class of digital filters that operate without 
arithmetic roundoff. These filters are linear, and can be used both as 
oscillators for signal generation and also as receivers for signal detec- 
tion. The feedback loop of each filter is constructed in such a way as 
to eliminate the possibility of roundoff or truncation errors, thus insur- 
ing perfect arithmetic. This entirely eliminates the problem of limit 
cycles. The filters presented, when used as generators, produce quan- 
tized values of A sin at of arbitrary accuracy. Implementation of these 
filters as receivers involves first sampling an analog input signal to 
produce a digital input into the filter. The filter is designed to resonate 
for a particular input frequency, thus enabling detection. 

The means by which arithmetic errors are eliminated in the feedback 
loop involves constraining all feedback coefficients to be integers (a 
constraint which turns out to be necessary to guarantee perfect 
arithmetic in any digital filter). Thus multiplication by these coeffi- 
cients can be performed as additions, simplifying implementation. 

The behavior of the feedback loop of this filter is modeled by a linear 
recursion whose characteristic polynomial is a cyclotomic polynomial. 
In recognition of this, we call the filter consisting of the feedback loop 
alone a "cyclotomic filter." It will be demonstrated that the only way 
to ensure perfect arithmetic with no limit on the period of operation 
(and thus avoid limit cycles) in a filter modeled by a linear recursion 
(i.e., a linear digital filter) is to constrain the feedback coefficients to 
be integers. Furthermore, it will be shown that, with this constraint 
on the feedback coefficients and also subject to minimizing memory 
and eliminating as many resonant harmonics as possible, the cyclotomic 
filter is uniquely optimal among all digital linear filters, both for the 
purpose of tone generation and the purpose of tone detection. 
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In subsequent sections of this paper, it is demonstrated how a 
weighting function can be applied externally to the cyclotomic filter 
to drastically reduce the impact of those higher-order resonances that 
remain. This is applied also to determine those impulse responses 
which have a small number of integer levels and lack higher-order 
harmonics. All the cyclotomic filters of practical significance, along 
with their associated weighting functions and impulse responses, are 
examined. 

In Ref. 3, a specific proposal is described for the Touch-Tone receiver 
(and tone generator), utilizing eight cyclotomic filters. 



II. CYCLOTOMIC FILTERS 

The purpose of this filter, as discussed in Section I, is to generate or 
detect a single pure tone u(t) = A sin (2irft + <p) of frequency /. 
Digital implementation involves realizing a discrete time filter with 
A; stages of memory (see Fig. 1), which is described recursively in terms 
of an input sequence u„ as 



= Y, aiX„-i + u n . 



(1) 



The numbers a»(* = 1, • • •, fc) are the feedback coefficients of the 
filter. The filter is driven by a clock with the time interval r between 
pulses. In tone generation, the filter must satisfy 



= u(nr), 



(2) 



at least for some initial conditions x , ■ • • , x k -i. When used as a 
receiver, the analog input u(l) is sampled, producing a discrete input 
Un = Asm {2rjiiT + <p) ; the filter (1) must distinguish between the 
desired frequency /o and all other frequencies in a band containing / . 
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Fig. 1 — Recursive filter in k stages of memory. 
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Specifically, it must satisfy the resonance property 

lim sup |.t„| = °o, (3) 

when / = /o, and in a sufficiently large band B including / there must 
be no other such resonances. Then \x n \ will be uniformly bounded in 
B in the complement of any small interval 5 about / , say, \x n \ ^ m(8) 
for all / E B, f ££ 8, for all n. A threshold detector can thus detect in 
a finite amount of time Nt, the presence (or absence) within B of an 
input frequency f (with error ±£|5|). It does this by comparing the 
gain sup „ $n \x n \ with the bound m (5) ; if sup \x n \ > w(5),then/G 5; 
otherwise it is not. Of course, the smaller the allowable error 5, the 
larger N must be. 

To know precisely when an input u n will resonate with respect to 
this filter, we first observe that the general solution to (1) is 

x n = t £ btf-%, (4) 

y-o <=i 

where pi, • • ■ , Pk are the roots (assumed to be distinct) of the charac- 
teristic equation 

X* - E a<x*-* - (5) 

and b x , • • • , b k are complex functions of the roots. [This is derived in 
(17) below.] If the magnitude of a root of (5) is greater than 1, the 
filter will be unstable. However, if all roots are inside the unit circle, 
then (1) will not have any resonance as defined in (3). Hence, in 
general we will assume that all roots of (5) lie inside or on the unit circle. 
Hence, the resonance (3) will occur if and only if the frequency / is 
such that with 6(i) = arg p, either 

2t/t - 6(i) (modulo 2tt) or 2tt/t = -6(i) (modulo 2») (6) 

for some i = 1, ■ ■ •, k with the property that \p%\ = 1. That is, the 
detector (see Fig. 2) will give a "yes" response iff (6) is satisfied. As we 
are trying to detect the presence of the frequency / = / , let us sup- 
pose by way of example that 6(1) = 2irf T (|pi| = 1). Then an input 
A sin (2irf t + <p) would elicit a "yes" response from our receiver. (Any 
phase shift of A sin 27r/ i will not affect the resonance of this signal, as 
A sin (2ir/ < + <p) = (A cos <p) sin 2irfot + (A sin <p) cos 2irfot, and 
cos <p and sin <p never simultaneously vanish.) However, let us now 
suppose that also 6(2) = 2h-/it (|p 2 | = 1). Then the receiver would 
also detect an input frequency /i (and would not differentiate between 
/o and /i). Hence, one would know only whether or not either / or /i 
is among the inputs. To positively identify the presence of / , one 
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Fig. 2 — Structure of a tone detector. 



must either insure that /i is out of band or use some other means to 
differentiate between /o and f\. 

Similarly, because of (6), the filter cannot distinguish between the 
frequency / and the frequency t -1 — /, since 2ir(r _1 — f)r = 2ir 
— 2tt/t = — 2t/t (modulo 2r) . In fact, 2irfr and — 2irfr are the re- 
spective arguments of complex conjugates, and thus we see from (6) 
that no new resonances can occur if the characteristic polynomial (5) 
is altered to include among its roots any complex conjugates of p\, • • • , 
Pk- We shall use this fact in our determination of a good structure for 
the recursion (4) . When the filter is such that an input of frequency 
/ will resonate, we shall say that the filter resonates (or has a reson- 
ance) at /. 

Recapitulating, because of (6), whenever the filter has a resonance 
at a frequency /, it will also necessarily and unavoidably resonate at 
the frequency t -1 — /. To counter the effect of this in practice, t must 
be made sufficiently small so that r -1 — / is out of band. In keeping 
with (6), we refer to resonance at the frequency / as "resonance at the 
root e i2r/T ," and resonance at t -1 — / as "resonance at the conjugate 
root e~ iUfT " [the roots in question being, of course, roots of (5)]. 

The remaining resonances described by (6) are those due to aliasing. 
These also are intrinsic to the system — a consequence of using discrete 
(rather than continuous) input samples u„. Indeed, if resonance occurs 
at a frequency / (or, equivalently, at the root e i2rfT ), it will also occur 
at all the frequencies / + mr~ l for any integer m, as 2rfr 
= 2ir(f + wit -1 )t (modulo 2x) or, equivalently, 

^/r = exp[i27r(/ + mf- l )r]. 

In practice, if conjugate resonances are out of band, resonances due to 
aliasing will also necessarily be out of band. 
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Hence, if pi, ■ • • , p m are those roots of (5) of modulus 1, the filter will 
have resonances within the band [0, t _1 ] at the frequencies 0(1)/2ttt, 
[2ir - 0(1)]/2ttt, 0(2)/2ttt, • • •, [2t - 0(w)]/2itt. The number of 
distinct resonances in the interval [0, (2t) -1 ] is m, less the number of 
roots among p\, ■ • ■ , p m which appear along with their conjugates. This 
picture is repeated in each successive interval [nr 1 , (n + 1)t _1 ] 
(n= ±1, ±2, ■•■)• 

It should be clear that, in choosing the recursion (1), one desires to 
have the number of resonances as small as possible — for the purpose 
of generation, to minimize the number of harmonics that can be pro- 
duced by perturbations of the initial conditions, and for the purpose of 
detection, to maximize the band in which the filter can detect a unique 
signal. Also, of course, one desires to have the memory k (a measure of 
the complexity of implementation) as small as possible. 

Ideally, one would like to have only one resonance, namely at the 
frequency one is trying to detect or generate. This is possible within 
the band [0, t -1 ], by using the recursion x n = —x n -i + u„. However, 
this resonates at a frequency equal to half the clock frequency t _1 and 
thus also resonates at the third harmonic (2t) _1 + t _1 due to aliasing. 
As the third harmonic is frequently in band, this recursion is generally 
not satisfactory. 

On the other hand, for some complex number p of unit modulus, one 
could use the recursion x n = px n -i + w„ which also has memory one. 
By adjusting t, one could make the argument of p = exp (i2irf T) 
small, thus avoiding any resonance up to as high a frequency as 
desired. However, there are problems with this recursion. First of all, 
the memory (in implementation) is not really one but two, as the real 
and imaginary parts of p must be handled separately. In fact, as seen 
before, no new resonances would be introduced by including the com- 
plex conjugate p of p to form a recursion of order two. Hence, one does 
just as well by replacing the characteristic equation X — p = with 
= (X — p) (X — p) = X 2 — aX + 1 (where the real number a = p 
+ p). The corresponding recursion replacing x„ = px n -i + u n , also 
(but now explicitly) of memory two, is x n = ax n -x — x n -2 + w„. This 
is the recursion after which digital linear filters are customarily 
modeled. However, as a (p) is, in general, not a rational number 
(gaussian rational*), it must in general be truncated, leading to slight 
frequency shifts, and multiplication round-off error in the feedback 
loop of these filters (Fig. 1) ; this could lead to unwanted limit cycles. 2 
To avoid this, a (p) is restricted to be rational (gaussian rational). 
Even for rational numbers, however, truncation error would occur if 



* Has rational real and imaginary parts. 
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the number of bits necessary to represent the number x n exceeded the 
word length allowed. In Section V we show that this can be controlled 
only if a is an integer. 

Hence, in the case of the real recursion, we restrict a to be an integer, 
and the only possibilities are = 0, ±1, ±2. We have already ruled 
out a = — 2 (this gives the square of the characteristic equation of 
x n +i = —x n + u n ). If a = 2, this gives the square of the characteristic 
equation of x„+i = x„ + u„, which is even worse, as it produces 
resonance at the second harmonic. The remaining three possibilities 
for a correspond to cyclotomic polynomials of orders 3, 4, and 6 (as 
defined subsequently in this section). It will be shown that, by taking 
a cyclotomic polynomial for the characteristic equation (5), one 
always obtains the best possible recursion (1) for the given amount of 
memory. 

In general, to have perfect arithmetic (the only means by which to 
uniformly avoid unwanted limit cycles), it is necessary to constrain the 
feedback coefficients a», i = 1, • • •, k [see (1)] to be gaussian integers 
(see Section V) . In fact, it will be shown that one can take each a* = 0, 
±1 so that each tap in the feedback loop involves at most changing 
the sign. Hence, from here on we restrict ourselves to three cases : the 
a/s are gaussian integers, are integers, or are 0, ±1. In what follows, 
we will show that the three are, for practical purposes, equivalent. 

For the first case, in the recursion corresponding to X — p = (no 
complex conjugate), the restriction to integer real and imaginary parts 
requires p = ±1, zki resulting in less generality than possible, as this 
corresponds to the recursions of the previous example with a — ±2 
only. In fact, we can generalize this, and say it is always better to 
include among the roots of (5) all the complex conjugates, and thus to 
have a recursion (1), all of whose coefficients are real (and hence 
integers). We will make this explicit in a moment, but let us first 
indicate the reasoning. First of all, by including the conjugates, no new 
resonances are introduced (as has already been demonstrated). Second, 
if among the roots of (5) even one conjugate were missing, the coeffi- 
cients of (1) would not all be real. In this case, the real and imaginary 
parts of x n would have to be considered separately, and one would thus 
need an effective memory of 2k. On the other hand, if one multiplies 
(5) by factors of the form (X — p), one for each root p of (5) whose 
complex conjugate is not also a root of (5), then the resulting poly- 
nomial and the corresponding recursion will have real coefficients. The 
respective degree and memory will thus be raised to no more than 2k 
(the effective memory of the complex recursion). Furthermore, as will 
be shown in Theorem 1 below, the new polynomial (and recursion) 
obtained from multiplication by the factors (X — p) will also be 
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guaranteed to have integer coefficients. Thus, we will do at least as 
well (and, as we have seen above, even better) by restricting all the 
recursions (1) to have real (and hence integer) coefficients. 
Let us now make this explicit. Suppose one has the recursion 



k 

%n = 2^ °l}%n—i I U n , 
3-1 



where aj(j = 1, • • •, k) are gaussian integers: ay = ay + b,i (ay and b } - 
integers, i = V— 1). Let y n and z„ be, respectively, the real and im- 
aginary parts of x n - Then 

k 

k 

3-1 

The only feature possibly mitigating in favor of the complex recursion 
is this: We are constrained to have ay and by be integers. If the new 
recursion with added roots did not have integer coefficients, then in 
spite of the other considerations above, one would choose the complex 
recursion. However, in the following theorem we show this is not 
possible. 

Theorem 1 : Suppose F(X) is a polynomial with gaussian integer coeffi- 
cients, and suppose pi, • • • , p m are those roots of F(\) whose complex con- 
jugates are not also roots of F. Then F(\) nf=i (^ ~~ Pi) has integer 
coefficients. Furthermore, if F(\) has no polynomial with integer coeffi- 
cients as a factor, then deg F = m. 

Proof: Write F(\) = g(\)h(\), where h(\) = J[ (X - p.-). Then g has 
real coefficients. Let p(X) be any irreducible factor of F(\) (considered 
as a polynomial over the gaussian integers). Suppose p has the root r 
in common with g and the root s in common with h. Then p (the 
polynomial in X whose coefficients are the complex conjugates of the 
coefficients of p) has f as a root, and hence p must also be a factor of 
F. But s is also a root of p, whereas s is expressly not a root of F. 
Hence, any irreducible factor of F must be a factor of either g or h. It 
follows that g has integer coefficients, and h (and thus h) have gaussian 
integer coefficients. As h(\)h(X) has real, and hence integer, coefficients 
the theorem follows. 

Thus, it is best to take the coefficients of the recursion (1) to be 
integers. The theorem which follows completely characterizes those 
recursions. 
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First, however, a short description of cyclotomic polynomials must 
be given. The Euler ^-function is a function on the positive integers, 
denned as follows: <p(m) is the number of positive integers less than or 
equal to m and having no integer factor in common with m, other than 
1 (such integers are said to be relatively prime to m). For example, 
<p(l) = v?(2) = 1, v>(3) - v(4) = 2, vK9) = 6. The cyclotomic ("cir- 
cle-dividing") polynomial of order m, denoted F m (\), is that monic 
polynomial (coefficient of the term of highest degree is 1) with integer 
coefficients all of whose roots are primitive mth roots of unity (that 
is, r m = 1, and r n ^ 1 for < n < m). Over the integers, F m (\) is 
irreducible (not a nontrivial product of polynomials with integer coeffi- 
cients). 4 From the definition, one can explicitly determine that F m (X) 
= lid (X — exp [_2iri{d/m)~\), where the product is taken over all d, 
1 ^ d < m such that d and m are relatively prime. Thus, the degree 
of F m is ip(m). 

The next theorem shows that, whatever constraints there are on 
available memory and acceptable resonant harmonics, the characteris- 
tic polynomial of the optimal recursion will be a cyclotomic polynomial. 

Theorem 2 : Let F(\) = X k - £* =1 a,X*-«', where a k ^ 0, a t (i = 1, • • • , k) 
are integers. Suppose every root p of F(\) = satisfies \p\ ^ 1. Then F 
is a product of cyclotomic polynomials. 

This is proved in Section V. Recall from our prior discussion that all 
the roots of F must be chosen to satisfy \p\ ^ 1 to have stable detec- 
tion. As it is, of course, better to have fewer resonances, one would 
hence choose for (4) a single cyclotomic polynomial. The cyclotomic 
polynomials make very desirable characteristic polynomials because of 
their extremely simple structure. For example, for m < 105 or for m a 
product of two primes, the coefficients of F m are all 0, ±1! For m a 
power of a single prime, the coefficients are all 0, 1 and for m < 385, 
the coefficients do not exceed 2 in absolute value. If m is a product of 
three distinct odd primes, all the coefficients are less than the smallest 
of those primes. These assertions are cited in Ref. 5. 

This means that implementation of the recursion (1) in the filter 
shown in Fig. 1 is very simple indeed. For all cases of practical interest, 
the feedback coefficients a, will be 0, ±1. Of course, when a, = 0, one 
simply does not put a tap on the ith stage. Because of the relation 

F pr ... P -(X) = F Pl ... Pn (X^r i --Kn 

(pi distinct primes — see Ref. 4), most of the coefficients of F m will 
usually be zero, and hence the taps-to-memory ratio is generally low 
(see Table I). 

In the preceding discussion, the principal emphasis has been on the 
use of the filter as a receiver. However, considerations relating to its 
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use as a generator lead to the same conclusion : that the characteristic 
polynomial (5) of the recursion (1) should be a cyclotomic polynomial. 
Indeed, for a generator, the problem of unwanted limit cycles is more 
critical. There is again the requirement that all the roots pi of (5) 
satisfy |p,| £» 1, as small perturbations in the initial conditions 
Xo, • • • , Xk-i from the (ideal) values 0, sin 2irf Q T, ■ • ■ , sin 2wf (k — l)r 
(to generate sin 2t f n t) are inevitable ; if such a perturbation occurs 
along an eigenvector corresponding to a root p„ where \pt\ > 1, it 
produces a nonzero coefficient 6, for that root in the general solution 
x n = £J=i bipl (where b h ■ • • , b k are functions of the initial conditions 
x 0) ••', x k -i; see Section III). This component would attain an arbi- 
trarily large amplitude (with time) and overwhelm the desired tone. 

Hence, one again requires a filter that can perform perfect arithmetic 
and whose characteristic equation has all its roots on the unit disc. 
From Theorem 2 we thus deduce that (5) should be a product of cyclo- 
tomic polynomials for the generator as well. As tone generation is 
impeded by the presence of harmonic resonances at other roots (due, 
again, to perturbation of initial conditions), one takes for (5) a single 
cyclotomic polynomial. 

Thus we have shown that, for both generating and receiving, the 
best linear recursion is one whose characteristic polynomial is cyclo- 
tomic. As the roots in this case are all of the form exp \_2iri(d/m)~\, the 
resonant frequencies can be expressed as 

2it/t = 2t — (modulo 2*) (7) 

m 

for all positive integers d < m such that d is relatively prime to m. 
Resonance at the fundamental is described by 2irfr = 2ir(l/m), that 
is, the fundamental of the filter is / = T~ l /m. Hence, if one requires a 
fundamental frequency of / (i.e., if / is the frequency of the tone to 
be generated or detected) and one intends to use a filter with memory 
k = <p(m), the clock rate r -1 is set at t -1 = f m. All other resonances 
occur at various harmonics (multiples of / ) as follows : the resonant 
harmonics in the band ^ / ^ t _1 occur when fr = d/m, that is, at 
f = df for all those integers d as above. For example, if m = 30 then 
k = 8 and d assumes the values 1, 7, 11, 13, 17, 19, 23, 29. Hence, this 
filter has no resonances between the fundamental / and the seventh 
harmonic. It resonates at the seventh harmonic 7/ , and thereafter at 
ll/o, 13/o, and so on. The resonances are at all the prime harmonics 
greater than 5, since in general those integers less than and relatively 
prime to the product m of the first p primes, are those primes lying 
between the pth prime and m. Furthermore, note that 30 = 1 + 29 
= 7 + 23 = 11 + 19 = 13 + 17. The first resonance due to aliasing 
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will always be at / = /o + t _1 = f + fom = (m + l)/ . In the case 
of the previous example, this is the thirty-first harmonic. 

Factors pertinent to the choice of which cyclotomic polynomial to 
use are relegated to Section VI. Suffice it to say at this point that the 
more memory available, the farther away from the fundamental can 
the first resonance be made due to aliasing. However, except for the 
cases m = 1 and m = 2, the first resonance after the fundamental will 
be below the clock frequency r~ l . In these cases, for a given amount of 
memory k, if the interest is to have the first higher-order resonance as 
far from the fundamental as possible, one would find the largest integer 
r such that the product m of the first r primes satisfies <p(m) ^ &• Then 
the first higher-order resonance would occur at the gth harmonic, where 
q is the (r + l)st prime. 

III. ELIMINATING IN-BAND HIGHER-ORDER RESONANCES 

The preceding analysis has indicated that, within the constraints 
established, various higher-order resonances are unavoidable. This 
could lead to difficulties. In practice, many higher-order harmonics are 
introduced in the process of limiting the input signal. The limiter (see 
Fig. 2) limits the amplitude of the input signal u(t). For example, a 
common limiter is a "hard-clipper." This has output ±1, depending 
upon whether u(t) ^ or u(t) < 0. The effect of hard-clipping on an 
input is to produce all the odd harmonics : sin 2vjt — > 2/ir sin 2irft 
-\- 2/37T sin Qirft + 2/5tt sin lQirft + • • • . Hence, a filter with more 
resonances frequently must be run for a longer period of time to attain 
a threshold sufficiently high to reject spurious signals. Also, when used 
as a generator, perturbations of the initial conditions of the filter could 
lead to unwanted harmonics at all the resonances of the filter. As such 
perturbations are inevitable, it is usually necessary to make allowance 
for eliminating these harmonics. 

While resonances due to aliasing are inherent to the discrete-time 
nature of the system and are hence unavoidable, resonances below the 
clock frequency t -1 can be handled outside the feedback loop. In par- 
ticular, it is possible (in theory) to eliminate (in practice, to reduce the 
Fourier coefficients of) any or all resonances at a frequency /, < / 
< (2t) _i , along with the conjugate resonance at t -1 — /. This is 
effected through operations outside the feedback loop. Specifically, this 
is accomplished either through alteration of the input before it enters 
the filter: u n — * v n = YJt=\ CiU n -i, or equivalently through alteration 
of the filter output before it enters the threshold detector: x„—>y n 
= Hi=i CiXn-i (see Figs. 3 and 8). Although these two options are 
mathematically equivalent, considerations with respect to minimizing 
the word length necessary for perfect arithmetic would mitigate in 
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favor of one or the other. This will be discussed in Section V. Here, we 
will describe the latter option only. 

Let X(X) be the generating function for the sequence 

k 

Xn = 2l< Q-iXn—i T W n , 
t"=l 

and let U(\) be the generating function for the input w„. That is, 

X(\) = t *n\ n , U(\) = £ u n \\ (8) 

n-0 n-0 

Then 

X(\) = £ o<X<X(X) + U(X), or X(X) = \ t/(X). 

'- 1 1 - £ a.V 

t = i 

Notice that denning F(X) = X* - £ a,X*-*', the characteristic poly- 
nomial of the filter, we obtain 

Since F(X) is assumed to be a cyclotomic polynomial, it is real and all 
its roots are of unit modulus. Hence p is a root if and only if p = p~ l is 
a root. It follows that X*F(X-') = F(X). Thus (9) may be rewritten as 

X(\) = ^ U(\). (10) 

We define a weighting function W (X) with the property that the 
resulting output function 

Y(\) m W(\)X(X) (11) 

has poles only at those roots of F(X) corresponding to those resonances 
actually desired. Specifically, W(X) will be a real polynomial of degree 
k — 2r, where r is the number of resonances desired in the band 
[0, (2r) -1 ]; the roots of W shall be those roots of F corresponding to 
the unwanted resonances. Typically, one desires to eliminate all 
resonances but the fundamental, in which case r = 1 and W(\)/F(\) 
= 1/(X 2 — aX + 1) for an appropriate real number a. Then, from (10) 
and (11), one obtains Y(\) = W(\)X(X) = (1/(X 2 - aX + l))t/(X) so 
Y(\) = -X 2 F(X) + a\Y(X) + U(\), and 

y n = ay„-i — y n -2 + u n . (12) 

This corresponds to a second-order filter with only one resonance in 
the band [0, (2t)- j ] as shown in Fig. 3. Although there will be trunca- 
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Fig. 3 — Implementation of the weighting function. 

tion error in (12), this will not lead to limit cycles, as there is no feed- 
back from this to the filter [although (12) represents the performance 
of the filter in terms of resonances, the filter, of course, is not realized 
in this way]. Specifically, the weighting function is implemented as in 
Fig. 3. This is derived from definition (11): if W(\) = £?=o c,X\ 
then equating terms in (11) yields 



y n = L CiXn-i, 

i=0 



(13) 



where, typically, d = k — 2. 

As mentioned earlier, the arithmetic of the weighting function is only 
approximate ; since there is truncation error in the computation of the 
coefficients c„ the roots of W will not precisely cancel out the roots of 
F. Rather, the roots of W will be slightly perturbed from the corre- 
sponding roots of F. The effect of this, as will be shown, is that all the 
resonances due to the roots of F (i.e., all the resonant harmonics of the 
original feedback loop) will be present in the output y n — however, they 
will have reduced energy (but for the fundamental). That is, the less 
the error in the implementation of W, the smaller the Fourier coeffi- 
cients of the higher resonant harmonics of the filter. This is demon- 
strated below. 

Suppose F is the cyclotomic polynomial of order m (or any poly- 
nomial whose roots p,, • • • , p* are distinct ???th roots of unity, so that 
each pj = e i2 * glm for some integer q, ^ q < m). A continuous-time 



GENERATOR-DETECTORS 481 



extension x(t) of the discrete-time function x„, satisfying x(ut) = x n 
can be denned as 

x(t) = E x n v{t - nr), (14) 

71=0 

where v(t) describes a continuous- time extension of x„. Specifically, 
v(t) is a periodic input pulse satisfying v(t + mr) = v(t) for all t 
[typically, v(t) = 1 for ^ t < r]. In (4), set w„ = v(wt) and nor- 
malize w(0) = 1. Then x n = E5=i M for n < m. Let f (q) [0(g)] 
denote the gth Fourier coefficient of x(t) [*>(0]- It follows that 

- v(q) E ^exp (-t'27r^w) 
n=o \ m / 

= Hq) E h E 1 P" exp (-*2* £ n) 

y=i n=o \ m / 

= Hq)bj, (15) 

where j is that index such that p, = exp [i2v(q/m)'}t if no such index 
exists, then x(q) = 0. To simplify matters, we will use the expression 
"the Fourier coefficient at (the root) p/' to indicate what in the case of 
(15) is the gth Fourier coefficient x (q). 

These Fourier coefficients can be computed explicitly from (9). 
Indeed, factoring X fc ^(\- 1 ) = II5=i (1 - pA) obtains 

jm\ 1 — pjh 

y=l 1 — pjh 

where the B/s are the coefficients of the partial fraction decomposition, 
derived explicitly in Lemma 3 below (it is assumed that all the roots 
Pi are distinct; in the case of multiple roots, however, similar results 
obtain). From (16) one obtains 

X(\) = EBit ( Pi X)« f Uik* 

i=l n— t'=0 

= E Bi E Pr'u**, (17) 

i *'.« 

so x n = E*=i B i E?=o p"~'wi [which is (4) above]. Hence, 5, = bj 
(j = i ( • •-,&) and their explicit form is given in the following lemma. 
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Lemma 3: Suppose pi, • • ■ , Pk are distinct numbers. Then 

A 1 ! 1 - Pi x kx V/=i (p. - Pi) ) i - p.-x 

;Vi 

Proof: The residue of the left-hand side at the ith pole is the coeffi- 
cient of that term in the sum above. The decomposition follows from 
the Cauchy residue theorem. 

Notice that, as the roots of F occur in conjugate pairs, a direct con- 
sequence of (17) is that, if p, and py are conjugate roots, then the cor- 
responding Fourier coefficients are also conjugate: &,• = fey. 

The Fourier coefficients for the sequence y n can be determined as in 
(15). For x„ = £ b/P* as before, we obtain from (13) 

y n = tcit btfT* 
i=o y=i 

= E Wfa)b jP % (18) 

y-i 

Thus, the Fourier coefficient of the sequence y„ at the root py is W(py)by 
(as could be expected, since Fourier transformations are multiplica- 
tive). Again, the conjugate coefficient W{p 3 ) = W(p 3 ). Observe that, 
if py is a root of W, then the Fourier coefficients of y n vanish at the roots 
py and py (W was chosen to be real). If W is the result of perturbing the 
coefficients of W to correspond to truncation error, then W'(py) is (by 
continuity) close to zero. Hence, as errors in the weighting functions 
are reduced, so is the power at each of the resonant harmonics above 
the fundamental (running the system for finite time, of course). 
Surprisingly, W is very stable ; if the coefficients of W are simply those 
of W rounded to the nearest integer (!), the results are frequently 
virtually as good as if W itself were used. This is exhibited in Table I 
and illustrated in Figs. 4, 5, and 6. These figures correspond to a filter 
using the cyclotomic polynomial F 30 . The input is a hard-clipped sine 
wave for each given frequency up to 15 times the fundamental. The 
input frequencies are normalized to units of the fundamental frequency 
for each filter. For each input frequency, the filter is run for an amount 
of time equal to seven cycles of the fundamental. If this time corre- 
sponds to Af steps of the filter, the output is max„ ^n \x n \, as measured 
at each input frequency (1500 samples). Using a W with integer 
coefficients (or any W with uniformly truncated coefficients) enables 
one to perform all the multiplications as additions, simplifying im- 
plementation and eliminating any further errors. As one expects, upon 
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FREQUENCY 

Fig. 4 — Hard-clipped/no weighting. 




4 6 8 10 12 

FREQUENCY 

Fig. 5 — Hard-clipped/rounded weighting. 



484 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1976 




FREQUENCY 

Fig. 6 — Hard-clipped/exact weighting. 



setting pi = pi, the Fourier coefficient of y n at the fundamental 



wm n 



pi 



pi 



A Pi(l " P;Pi) 



Pi 



y-a Pi — P> Pi — P2y=3 Pi _ Pj Pi _ P2 

is the Fourier coefficient of (12) at the fundamental. 

IV. IMPULSE RESPONSE 

The impulse response is the output resulting from an input of a 
single pulse: u = 1, u n>0 = 0. Since this output can also be produced 
by appropriately setting initial conditions, we will refer to it as a pulse 
train. From (4) we see that if the input u n is a single pulse, then the 
output x n reduces to 



= L b iP l 

t=i 



(19) 



In the context of the previous sections, it is assumed that the charac- 
teristic polynomial of the sequence x n is cyclotomic. Since each p.- is 
then an wth root of unity, the sequence x„ is periodic : x„ +m = x n for 
all n. As before, the resonant harmonics present in the pulse trains x n 
correspond to the mth roots of unity which are roots pi (i = 1, • • •, k) 
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of F = F m ; the Fourier coefficient of the pulse train at the root p, is 
bj (see Section III). 

In particular, using the notations of Section III, U(X) = 1 and thus 
(10) reduces to 

xm = ~ (20) 

But X(\) = T.n^0X n \ n = (E?-0 Xn\ n )CZn=0^ mn ) SUlCe X n+m = X„. 

Defining /(X) = X)»-o x n \ n , one obtains 

/W-^f (2D 

from (20). Notice that /has integer coefficients (the input u n is integer, 
as are the coefficients a,-). Indeed, 1 — X m is a product of cyclotomic 
polynomials, one of which is F(\). Specifically, 

1 - X- = ± II ^n(X) 
n\m 

[the product is taken over all n which divide m; hence, for example, 

1 _ x « = -F l (\)F 2 (\)F 3 (\)F 6 (\) 

= (X - 1)(X + 1)(X 2 + X + 1)(X 2 - X + 1)]; 

and, from (21), 

/(A) = ± n F n (\) 

n | m 

obtains. Consequently, /(p) = for all mth roots of unity p, except 
for the primitive roots of unity [the roots of F m (X)^}. This was antici- 
pated by E. N. Gilbert in Ref. 6, where he showed that a pulse train 
x„ of period m has resonances at those harmonics corresponding to the 
mth roots of unity which are not roots of YJn=l Xn\ n = 0. Equation 
(21) covers the general situation where /(X) [and consequently F(X)2 
are arbitrary products of cyclotomic factors of 1 — X m . 

In the same paper, Gilbert was concerned about the problem of 
increasing the power of the pulse train at the fundamental (relative to 
the power at the other resonances) . This could be done by shaping the 
input u n for one period, but it is usually undesirable to do this. As 
explained in Section III, however, the same effect is obtained by 
utilizing a weighting function W. If utilized directly, this will introduce 
noninteger levels into the pulse train. Nonetheless, it is possible to 
avoid this by replacing W with W 1 where the latter is obtained through 
rounding off to the nearest integer the coefficients of the former. The 
pulse train resulting from W 1 will have integer levels, but the trunca- 
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tion error will again introduce higher-order resonances. However, 
Table I shows that these are very small indeed, leaving typically about 
98 percent of the power at the fundamental. This compares with 25 
percent or less (for F u , F 2i , F 3 o) without W. Note, for example, F 8 . 
The pulse train 1, 0, 0, 0, — 1, 0, 0, has resonances at the third, fifth, 
and seventh harmonics. However, by simply altering this to 1, 1, 1, 0, 
— 1, —1, —1, 0, the first appreciable resonance does not come until 
the seventh harmonic. In this case, use of W 1 does not introduce any 
new levels in the pulse train. 

The worst case in Table I is F 9 where 92 percent of the power is at 
the fundamental. E. N. Gilbert has pointed out that if one wished to 
increase the proportion of the power at the fundamental of this train 
(or any other), one could multiply the output Y(\) by some constant 
c > 1, chosen so that the roundoff error of cW — > (cW) 1 is smaller 
than that for W 1 alone [recall (11)]. This, however, would introduce 
more levels into the pulse train (although no more than c times as 
many) . 

Table I gives an indication of the possibilities for various niters. 
Included are the niters with memory less than 12 which provide the 
greatest separation between the fundamental and the first resonant 
harmonic, either with or without the weighting function. The asterisks 
and daggers indicate those which, for the amount of memory, have the 
largest possible separation without or with the weighting function. For 
utilization with a "hard-clipper" (which has all odd harmonics), F 3 , 
Fg, and Fi 5 are included. Although these resonate at all even harmonics, 
they have the same response to a hard-clipped input at the fundamental 
as the respective cyclotomic filters of twice the sampling rate. To have 
the first resonant harmonic higher than the seventh (without W) 
would require a memory of 48 (and F to have a coefficient of 2). The 
next interesting entry with respect to W is F 36 with memory 12. The 
columns to the right of the double line all deal with the integer-rounded 
transfer function W 1 . Columns A and B give |bi| 2 /L*=i |6.-| 2 and 
|6iW / (pi) | 2 /L*=i \biW z (pi)\ 2 as a percent, respectively, where &,• is 
the Fourier coefficient of the sequence x„ at the root p* [see (6) and 
Section III)]. Column C gives (max 2 g,g* /2 | W(p.) | 2 )/| W(pi) | 2 as 
a percent. The roots p, . (i = 1, • • • , k/2) are assumed to be in order of 
ascending argument <ir (so pi is the fundamental). Columns D and E 
give the moduli of the Fourier coefficients 6, and b i W I (pi) of the 
sequences x„ and y„, respectively. Columns F and G give the pulse 
trains of x n and y„, respectively, with initial pulse uo = 1, u„ >0 = 0. 
The exponent denotes repeated digit; the arrow indicates that the 
preceding train is followed by another identical train, but that each 
digit is the negative of what it was. 
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Table I — Characteristics of 



Mem- 
ory 


Characteristic 
Polynomial 

F 


Resonant Harmonics in the 
Band [0, r -1 ], Aside From 

Fundamental 


First 

Resonant 

Harmonic Due 

to Aliasing 


Number of 

TapB on Filter 

(Without 

Weighting 

Function) 


Integer-Rounded 
Weighting 
Function 
W 




Without 
Weighting 
Function 


With 
Weighting 
Function 




*1 
2 
2 

*2 
4 

n 

6 

te 

8 

8 

8 

t*8 


F, = X+1 

Fi-X'+X+l 

F.=X»+1 

F«=X»-X+1 

F B =X<+1 

F«-X«-X«+l 

Fs=X«+X»+l 

Pu = X«-X'+l 

Fi5 = X»-X'+X 6 
-X<+X"-X 

+1 

Fi«=-X8+1 

Fi« = X»-X«+l 

Fio = X«+X»-\« 
-X«-X»+X 
+ 1 


none 

2 

3 

5 

3, 5, 7 

5, 7, 11 

2, 4, 5, 7, 8 

5, 7, 11, 13, 17 

2, 4, 7, 8, 11, 
13, 14 

3, 5, 7, 9, 11, 
13, 15 

5, 7, 11, 13, 17, 

19, 23 
7, 11, 13, 17, 

19, 23, 29 


7 
11 
8 

17 

14 

15 
23 

29 


3 
4 
5 
7 
9 
13 
10 

19 

16 

17 
25 
31 


1 
2 
1 
2 
1 
2 
2 

2 

6 

1 
2 
6 


1+X+X" 

1+2X+X' 

1+2X+X»+2X» 

+X« 
1 + 2X+3X»+2X" 

+x« 

1+X+X»+X«+X« 
+X»+X« 

l-r-2X+2X 1 +3X> 
+ 2X«+2X6+X« 

1+2X+3X»+3X» 
+ 3X«+2X«+X» 

1+3X+5X»+5X> 
+ 5X'+3X»+X« 



*t See text for explanation 

V. CONDITIONS FOR PERFECT ARITHMETIC 

Here we indicate why cyclotomic polynomials yield optimal recur- 
sions for generating sinuosidal signals. When we use (1) to generate 
tones, the u n is set to zero and some initial condition x , Xi, • • • , Xk-i is 
chosen to generate the required samples x„ : 



«-l 



(22) 



If we use the usual second-order recursion, then (22) is of the form 

X n = aXn-1 — Xn-2, (23) 

where \a\ < 2, so we have complex roots. In this case, we show below 
that the number of distinct values that x n , n = 0, 1, • ■ •, N can take 
grows at least as fast as N/2, with N. So, to simulate (22) with perfect 
arithmetic, the number of "words" needed grows at least as fast as N, 
the number of samples needed. 

Proposition 4 : Suppose \a\ < 2, and rational but not an integer. Then 
for any initial conditions x Q , xi (not both zero) and any positive integer N, 
the number of distinct values among xo, • • • , xn, where x n = ax n -\ 
— z„_2, for 2 ^ n ^ N, is at least N/2. 
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some 


cyclotomic filters 










A 


B 


Highest Power 


D 


E 


F 


G 


% of Total 










Power in 
CO, (2t)-'] at 
Fundamental 


at a Rejected 

Resonance 

(%of 


Coefficients (in Order of 
Arguments < x) 


Pulse Trains 




Fundamental) 










Without 
W 


With 
W 


Xn 


1/n 


Xn 


Vn 


100 








1 


— 


1, -1 





100 


— 


— 


0.58 


— 


1, -1,0 





100 





— 


0.5 


— 


1.0, -1,0 





100 


— 


— 


1.73 


— 


1*0—1*0 





50 


97.1 


2.9 


All 0.25 


0.60, 0.10 


10'- 10' 


1»0-1>0 


50 


99.5 


0.5 


All 0.29 


1.08, 0.08 


1010' -10 -10' 


12«10-l-2'-10 


33.3 


92.7 


7.8 


All 0.19 


0.85, 0.04, 
0.24 


10*- 10* 


121*-l*-2-10 


33.3 


99.6 


0.3 


All 0.19 


7.6, 0.04, 
0.06 


iono s — 


123*210- 


53.3 


98.9 


0.8 


0.33, 0.27, 
0.09, 0.11 


4.78, 0.51, 
0.55, 0.75 


1«0»— 1«0» 


123*21-1-2 
-3«-2-10 


25 


98.0 


1.8 


All 0.13 


1.29, 0.02, 
0.06, 0.17 


10'- 10' 


12*32*10- 


25 


99.5 


0.3 


All 0.14 


1.97, 0.05, 
0.09, 0.11 


10'10'-10'-10' 


123*4>3»210- 


6 


98.4 


1.0 


0.11, 0.09, 
0.27, 0.33 


2.41, 0.04, 
0.19, 0.24 


i - lion- no' -♦ 


123»545*453»210- 



Proof: We can write x n = bip" + &2P2, where pi, pi are the distinct 
roots of X 2 — a\ 4- 1, as in (19). Since the roots are not real, let 
P = pi(= Pi), b = bi(= 6 2 ). Then z„ = x m implies Re (bp n ) = Re (6p m ). 
In this case, letting = arg p, tp = arg b, we obtain cos (<p + nd) 
= cos (<p + md) so ip + nd = ±(^ + md) (mod 2r). Since p is not a 
root of unity, the numbers nd (n = 0, 1, 2, • • •) are all distinct and 
hence for fixed in either n = m or nd = —2<p — md (mod 2t). As this 
last congruence can be satisfied by at most one n, it follows that, for 
each ;/(, there is at most one n 5* m such that x n = x m . 

The following result shows that, if one wishes to generate sin rnB 
with perfect accuracy using a linear recursion, e irB must be a root of 
the corresponding polynomial (5). 

Proposition 5: If s n = sin irnd is a solution of x„ = £ ajX n -j and 6 is 
not an integer, then e" e is a root of the polynomial \ k — £ aX* -3 . ' 

Proof: From sin ir{n + 1)0 = L ay sin ir(n + 1 - j)d, we expand 
both sides using a familiar trigonometric identity and get 

sin irnd cos ird -f- cos trJid sin ttd = sin ir(n + 1)0 

= 2 a i sin ir(n + 1 — j)d = 5Z <*j sin ir(n — j)d cos ir0 
+ 2 aj cos ir(n — j)6 sin 7r0 = sin irnd cos ird 

+ 2Z Q-i cos ir{n — j)d sin ird. 
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Since d is not an integer, sin ird ^ 0, and thus from the equality of the 
first and last expressions, we obtain cos trnd = £ <*/ cos v(n — j)d. 
Hence, cos irnd is also a solution to the recursion, and it follows that 
giwnt = cos xne _|_ { s i n Tn o [ s a solution too. Consequently, e" M 

- E a^'X*-')* = 0. 

The next theorem shows that every recursion which satisfies the 
stability criterion \p\ ^ 1 for all its roots, and for which perfect 
arithmetic is possible, is cyclotomic. 

Theorem 6: Suppose every root p of the polynomial F(\) = X* 

— £*=i a,-A* - ' satisfies \p\ ^ 1. 

(i) If ai, • • •, a k are integers and a k 5* 0, then F(\) is a product of 

cyclotomic polynomials, 
(ii) If Oi, • • • , a k are rational numbers and x n = Y, a-i^n-i is periodic 
(Xn+p = x n for some p, all n) for some nonzero initial conditions 
x , •••, x k -i, then F(\) has as a factor a cyclotomic polynomial. 

Proof: For case (i), each irreducible factor (over the integers) of F(X) 
has the same form as F(\) itself by' 'Gauss' Lemma". 6 Thus, it 
suffices to assume that F(\) is irreducible, in which case all its roots are 
distinct. In this case, we can write x n = £ 6,-p? where the p.'s are the 
roots of F(\) and x„ is as in case (ii). But then \x n \ ^ L I °i I » and as 
for any integer initial conditions x , • • ■ , x k -i, x n will be an integer for 
all n, x n can in such a case assume only a finite number m of distinct 
values (m = [2ZI^«I!])- Hence for all n, the fctuple (.x n +i, • ■ •, x n +k) 
can assume at most m k distinct values, and as x n is recursively gen- 
erated with memory k, x n must be periodic, of period p ^ m*. This 
brings us to case {ii). 

For case (ii), let L be the rational canonical form associated with 
the recursion x n (see Ref. 7, Section 5.2.1), and J be the Jordan 
canonical form of L. Then for some initial state vector x, J"x = x, and 
it follows that some diagonal element of J, that is, some root of F(\), 
must be a pth root of unity. Hence, the irreducible factor of F(\) 
having that root must be cyclotomic. 

Hence, from the above the 6 of Proposition 5 must be rational when 
perfect accuracy is required. 

In all the preceding, the basic assumption has been that all the 
coefficients of the recursion (22) are real. We can infer from Theorem 1 
that this is no loss of generality as, if the recursion had complex coeffi- 
cients (with rational real and imaginary parts) and was irreducible 
over the field Q(i) (the field of gaussian rationals), then the roots of 
the characteristic polynomial would be distinct, no pair being con- 
jugate. Indeed, Theorem 1 remains true if the word "integer" is every- 
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where replaced by "rational number." The arguments of Section II 
show that we may as well assume all the coefficients are real. 

VI. COMPUTING WORD LENGTH AND ADDITIONS PER CYCLE 

To realize the cyclotomic filters in hardware with perfect arithmetic, 
the necessary amount of memory and adder complexity must be pro- 
vided. We describe here how to estimate the word length and the rate 
of additions required to implement a cyclotomic filter with a weighting 
function. It shall be assumed that all operations are performed in 
binary form. The number of binary bits required to store each x n is 
called the word length o> of the system. For generators that produce a 
signal approximating a sinusoid, the word length required will depend 
on the accuracy of approximation needed. When the filter is used as a 
tone detector, the word length required will depend on the duration of 
operation, since the signal level tends to build up, especially at fre- 
quencies close to any resonant frequency (Fig. 7). The signal level, of 
course, does not uniquely specify the minimum word length. Even 
though for storing x„ we may need only w bits, it is conceivable that 
during the computations numbers greater in magnitude than x„, which 
need more bits for storage, could arise. To perform operations in a 
serial-multiplexed fashion, it is desirable to have uniform word length 
for all operations in the feedback loop of the filter. Hence, the word 
length will have to be increased to accommodate any number en- 
countered during the computations. However, for the filters considered 




56 60 

STEPS 

Fig. 7 — Growth of output. 
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in Table II, it is possible to arrange the computations in such a way 
that the word length is determined by the maximum magnitude of x n . 
In general, there are a finite number of ways in which the additions 
involved in the filter can be arranged. By simulation of the different 
arrangements, the word length required can then be determined. 

There are two possible ways of implementing the cyclotomic filters 
as generators. The first is to generate the impulse response (19) ; this 
is generally sufficient (see Table I). In this case, the weighting function 
(13) shapes the effect of this impulse to simulate the initial conditions 
x , • • •, x k -\ of the tone being generated. As the input is zero after the 
initial pulse u = 1, the weighting function need only be used during 
the first d + 1 steps of the filter. Let m be the largest number in the 
pulse train y n of Table I, and let [[>]] be tne smallest integer larger 
than x. The word length necessary for perfect arithmetic is at least 
03 = CD°g2 w~D + 1 an d, f° r * ne filters considered here, a> is also 
sufficient. (We add 1 for a sign bit.) This word length is shown in 
column B of Table II. 

However, rounding off in the weighting function introduces errors 
in the effective initial values of the signal. If this approximation is not 
sufficiently good, then the initial conditions of the filter x , • • • , Xk-i 
can be set as accurately as needed, and then the filter is operated with 
the feedback loop alone. In particular, one can set the initial conditions 
of the filter such that \x n - sin 2irn/p\ < 2~ m (n = 0, ■ ■ ■ , k - 1) 
where sin 2-irn/p is the desired signal. One can then compute the mini- 
mum word length required by simulating the filter for one period. In 
all cases of interest here, the word length including sign is (m + 1) 
for m ^ 12. Hence, as an example, the cyclotomic filter of order 30 
can generate a sequence (x„) such that \x n — sin 2irn/p\ < 2 -10 if 
the initial conditions are set such that \x n — sin 2irn/p\< 2 -10 
(n = 0, ■ • •, 7), using a word length of 11. 

To determine the number of binary additions per period of the filter 
(i.e., per cycle of the fundamental), one counts the number of bit 
additions per step. If m denotes the number of additions per step, then 
pmw is the number of binary additions per cycle, where p is the period 
of (22) and a> the word length used in the feedback loop (see above). 
When the generator is implemented in the first way (using an initial 
pulse and the weighting function), the number of additions is shown 
in column C of Table II (not including those necessary in the initial 
d + 1 steps for the weighting function). When the generator is im- 
plemented in the second way (setting the initial conditions), the num- 
ber of additions can be computed by multiplying the value in column C 
by w/u', where co is the word length chosen and co' is the corresponding 
word length from column B. 
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When the filter is used as a detector, we assume that the input to the 
filter is a sequence which only assumes the values +1 and —1. This is 
true, for example, when the analog signal to be detected is either hard- 
clipped or delta-modulated. In these cases, it is advantageous to apply 
the weight function to the input sequence u n rather than to the se- 
quence x n ; since, in general, x n can assume many values other than +1 
and —1, computations involving the weighting function are simplified 
if they are performed on the input (see Section III). In fact, applying 
the weighting function to the input is so simple arithmetically that it 
can be implemented with read-only memory. On the other hand, if 
read-only memory is not used and one wishes to save on computations 
by checking the threshold (max |x»J) only in the last cycle of the filter 
(with respect to its duration of operation for detection), then the 
weighting function is best implemented as in Section III, on the output 
of the feedback loop. Then the filter can be run during all but the last 
cycle, without computing the weighting function. 

When the weighting function is applied to the input, the filter is 
described by 



Vn = H C,W„_i 
t-Q 



(24) 






(25) 



where u„ is the input into the filter and v n is the result of the weighting 
function. Figure 8 describes this filter. 

For the filters in Table I, the effect of rounding c,- to the nearest 
integer is slight. Hence, it is a fortiori suitable to round off 
v n = L CiUn-i to the nearest integer. Therefore, since the only values 
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Fig. 8 — Implementation of the weighting function at the input. 
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assumed by w,- are ±1, it suffices to have for v n a word length of 
o» = [[log 2 {Y, \ci\ }]] + 1 (where {a;} is the integer closest to x and 
[[#]] is the smallest integer larger than x ; 1 is added for a sign bit) . 
The sequence v„ can then assume any value between — { £ | c, | } and 
{ T, I Ct | } . With d as in (24) and w as above, implementations of the 
weighting function with read-only memory then requires 2 d+1 oi memory 
bits. The respective values for this are shown in column D of Table II. 
When a bank of such tuned niters is used in one receiver (for example, 
in a Touch-Tone® system such as described in Ref. 3), all the filters 
could use one read-only memory for the weighting functions. Also, by 
increasing oj, we can make the round-off error as small as we wish. 

To determine the word length for use in the feedback loop of the 
detector, the maximum signal level can be determined by using an 
input u n of the same frequency as the resonant frequency. Since the 
impulse response [see (19)] of these filters is periodic and of the same 
period as the resonant frequency, the latter produces the maximum 
signal level sup n ^Ar x n , for duration of operation Nt. Let this maximum 
be M. The word length required should then be at least [[log 2 MJ] + %• 
For all the filters considered here, [[log 2 M~\~] + 1 is also sufficient. 
The number of M, of course, is determined by N. If the cyclotomic 
filter is of period p (i.e., Theorem 1 is F p ), then the filter runs through 
N/p periods, corresponding to N/p cycles of the fundamental. Calcula- 
tions have been made for two values of N/p : 7 (the number of cycles 
computed in Ref. 3 to be necessary for Touch-Tone interchannel rejec- 
tion), and 10 (a more uniform point of reference). 

In Table II, column E shows the word length required in the feed- 
back loop for the indicated durations, when the weighting function is 
computed on the input as in (24), implemented equivalently with or 
without read-only memory, producing the filter response (25). 

When there is no weighting function on the input, the word length 
required is shown in column F (of course, a weighting function may be 
applied to the output as in Section III). 

The number of binary additions per cycle for the detector is de- 
termined in the same way as for the generator ; the number is pmu as 
defined above. These numbers are shown in columns G, H, and K of 
Table II. Column G shows the number of binary additions per cycle 
in the feedback loop when read-only memory is used to implement the 
weighting function, applied to the input as in (24). If read-only memory 
is not used, then the weighting function has to be computed. Since the 
numbers involved in the computation of the weighting function [when 
implemented as in (24)] are generally smaller than those in the feed- 
back loop, the word length required for their computations are smaller. 
Hence, one can use two different adders, one for the weighting function 
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and one for the feedback loop. Using this arrangement, the number of 
additions per cycle for calculating the weighting function is shown in 
column K. The number of binary additions per cycle when no weighting 
function is used is shown in column H. This, of course, applies when 
the weighting function is applied to the output as in Section III (but 
does not include the number of additions necessary for the weighting 
function). To calculate the number of additions when the weighting 
function is applied to the input, but read-only memory is not used, add 
columns H and K. 

Column A indicates the respective cyclotomic niters described by 
their periods. 

One important consideration that affects the choice of the order of 
cyclotomic filter is the noise level at the input to limiter (together with 
the noise in the limiter) . This affects the output of the limiter when the 
signal level is low. One could divide the period of the signal to be 
detected into regions where errors could affect the decision about the 
sign of the signal, and regions where no errors will occur. Those sam- 
pling instances where errors could occur lie in regions where the absolute 
value of the signal is small. Suppose these regions are intervals of 
length e around the zero crossings of the signal. The worst case cor- 
responds to a phase shift of the signal with respect to the sampling 
interval which maximizes the number of samples in the error regions. 
For e = 1/63 (corresponding to approximately 20 dB s/n), there are 
at most two samples per period that are subject to errors for all the 
filters we have considered here. Hence the ratio of error-susceptible 
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samples to error-free ones decreases in this case as the period p in- 
creases (for p ±5 30) . This ratio indicates the perturbation of the thresh- 
old one has to make in order to compensate for errors in the limiter. 

VII. APPLICATIONS 

Possible uses for the systems described in this paper have been 
mentioned in Section I. In particular, a scheme is proposed in Ref. 3 
for utilizing eight cyclotomic niters as channel detectors in a Touch- 
Tone receiver. 

Another application of cyclotomic filters may be fsk. As described 
earlier, by selecting the initial conditions of a cyclotomic filter of 
period p, one can approximate uniformly sampled values of a sinusoid 
of period p, i.e., sin 2irn/p. By changing the clock rate of the filter, one 
can shift the frequency of the sinusoid to any preassigned value. 
Hence, when using the filter as a generator, one can shift the clock rate 
to shift the frequency. This method of shifting frequencies does not 
introduce any "discontinuities" in the signal. If, instead of changing 
clock rate, one were to change the coefficient of a filter, then the filter 
has to be reinitialized to have constant amplitude, thus producing a 
discontinuity in the signal. In a similar manner, when using the filter 
as a detector, one can shift the resonant frequency by shifting clock 
rate. Hence, with the same filter, one can generate and detect both 
tones used in a typical fsk arrangement. Furthermore, cyclotomic 
filters have infinite Q, allowing for the possibility of increasing signaling 
rate above the presently used systems with finite Q. 

REFERENCES 

1. L. B. Jackson, "An Analysis of Limit Cycles due to Multiplication Roundoff in 

Recursive Digital Filters," Proceedings of the Seventh Allerton Conference on 
Circuits and Systems Theory, 1969, pp. 69-78. 

2. H. Breece, private communication. 

3. B. Gopinath and R. P. Kurshan, "A Touch-Tone® Receiver-Generator with 

Digital Channel Filters," B.S.T.J., this issue, pp. 455-467. 

4. S. Lang, Algebra, New York: Addison-Wesley, 1965. 

5. E. Lehmer, "On the Magnitude of the Coefficients of the Cyclotomic Polynomial," 

Bull. Am. Math. Soc, 42 (June 1936), p. 389. 

6. E. N. Gilbert, "Pulse Trains Which Lack Prescribed Harmonics," unpublished 

document. 

7. B. Gopinath and R. P. Kurshan, "Recursively Generated Periodic Sequences," 

Canadian Journal of Math., XXVI, No. 6, 1974, pp. 1356-1371. 



496 THE BELL SYSTEM TECHNICAL JOURNAL, APRIL 1976 



